Relativistic mean field description 
for the shears band mechanism in 84 Rb 
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For the first time, the Relativistic Mean Field (RMF) model is applied to the shears band recently 
observed in j^Rb. Signals of the appearance of the shears mechanism, such as smooth decreases of the 
shears angle and of the B(M1)/ B(E2) ratio with keeping the nearly constant tilt angle, are well reproduced. 
Thus it is shown that the microscopic RMF model can nicely describe the shears band in this nucleus. 
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Recently obtained data of the so-called shears bands 
in the proton rich Pb-isotopes |l|-[5| are well described 
within the framework of the Tilted Axis Cranking (TAC) 
P-|TT| approach. In the shears bands, the magnetic dipole 
vector, which arises from few proton particles (holes) and 
few neutron holes (particles) in high-j orbitals, rotates 
around the total angular momentum vector. At the band 
head, the proton and neutron angular momenta are al- 
most perpendicular. This coupling results in the total an- 
gular momentum which is tilted from the principal axes. 
With an increase of the rotational frequency, both the 
proton and neutron angular momenta align toward the 
total angular momentum. Consequently, the direction of 
the total angular momentum does not change so much 
and regular rotational bands are formed in spite of the 
fact that the density distribution of the nucleus is almost 
spherical or weakly deformed. These kinds of rotation 
are called magnetic rotation in order to distinguish 

from the usual collective rotation in well-deformed nu- 
clei (called electric rotation). Magnetic rotation has also 
been observed in other regions such as A ~ 110 HHIl] 
and 140 p9| j regions. In a recent experiment, new shears 
bands in the Rb- isotopes were discovered 20 2^]. These 
are the first experimental data obtained in the A ~ 80 
mass region. From the theoretical side, such shears bands 
have been well examined by the shell model p3l and a 
model based on the mean field approximation^ . The 
shell model approach is especially suitable for those in 
the vicinity of doubly closed nuclei. When we add more 
and more particles to the doubly closed core, a transi- 
tion from complicated mult iplet- like level structures to 
regular rotational bands in deformed nuclei is observed. 
Shears bands are seen in the middle region, that is, the 
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system in which only few proton particles/holes and few 
neutron holes/particles are involved. The shell model 
can describe well both a multiplet-like and a rotational- 
like structure, including the shears bands, as well as a 
transition from one to another structure. 

A weak point of the shell model approach is the effects 
coming from the truncation of the model space. As for 
the study of shears bands, it is suggested that the 
polarization of the low-j orbitals (e.g. p/-shells in the Pb- 
isotopes) would have a 'glue' role, which keeps the outside 
high-j particles or holes being a blade of shears. This 
is important for the bands to be stabilized. The trun- 
cated shell model may not be suitable to describe such 
effects properly because of the limitation of the model 
space. On the other hand, such polarization effects can 
be easily included in the mean field models as there is 
no limitation of the configuration space. It is difficult, 
however, to describe the multiplet-like level structures 
and the transition from multiplet-like to rotational-like 
structures in the simple mean field approach. Therefore, 
complementary approaches, from both the shell model 
and the mean field models, are necessary in order to get 
a whole description of the property of the shears bands. 

Focusing on the mean field approaches to the shears 
bands, only studies based on the pairing+QQ Hamilto- 
nian |8|,|| have been done up to now. Those based on 
more sophisticated models, such as the Skyrme Hartree 
Fock (SHF) and the Relativistic Mean Field (RMF) mod- 
els are still missing and strongly desired. In a previous 
preliminary work [|23| , we have applied the RMF model to 
the tilted axis rotation. In the present paper, we for the 
first time describe the shears bands by the RMF model, 
which has gratefully been successful in reproducing many 
properties of finite nuclei. 

In the RMF model [gi] , we consider the following La- 
grangian, which contains the nucleon and several kinds of 
meson fields, such as a-, uj- and p-mesons, together with 
the photon fields (denoted by A) mediating the Coulomb 
interaction: 
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here £; n t is the interaction part between nucleons and 
mesons. The non-linear self interactions among the <j- 
mesons, £nlo-, are also included. 

For applications to rotating nuclei, the Lagrangian of 
the RMF model must be generalized into a uniformly 
rotating frame which rotates with a constant rotational 
frequency, CI = (Cl x ,Cl y ,Cl z ). From this generalized La- 
grangian, the equations of motion are derived. This can 
be done in the same manner as that in the Principal Axis 
Cranking (PAC) case 1 25 - 2^,|2j| . The resulting equations 
are 
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where the p-meson and photon fields are omitted for sim- 
plicity, although they are included in the numerical cal- 
culation. It is known J2?]] that the Coriolis terms for 
the meson and photon fields (those proportional to CI 2 
in the above equations) give very small contribution and 
can be completely neglected. The method used to solve 
the coupled equations of motion is, again, the same as 
that in the PAC case. The nucleon and meson fields are 
expanded in terms of 3-dimensional harmonic-oscillator 
eigenf unctions. Note that, contrary to the PAC case, 
the signature symmetry is broken in the TAC approach. 
The parity is thus the only symmetry we assume in our 
code. The cutoff parameters of the expansion are taken 
as Np — 10 for nucleon fields and Nb = 10 for meson 
fields, respectively. When we increase these cutoff pa- 
rameters to N-p = N-q = 12, we find the changes of the 
calculated values are only 0.1% for the total energies, 
2 — 3% for the total quadrupole moments and 0.5% for 
the total angular momenta, respectively. As for the pa- 
rameter set, we use that called NL3 |2{|. In the present 
code, the pairing correlations are not taken into account. 
They should be included for a precise description of the 
properties of heavy and medium-heavy nuclei in the low 
spin region. In the relativistic case, however, the size 
of the Hamiltonian in the Hartree-Bogoliubov equation 
becomes twice as large as that of the non-relativistic 
case due to the degrees of freedom of the lower com- 
ponents. This makes it very time-consuming to perform 
a 3-dimensional cranking calculations with pairing. Be- 
cause this is the first RMF work for the shears bands, 
one of the important purposes of which is to examine its 
applicability to the shears bands, we here concentrate on 
the calculation without including the pairing interaction. 

As the first example of our RMF calculation of the 
shears bands, we choose the nucleus ||Rb47. In all of 



our calculations, we assume that the proton configura- 
tion is fixed to be 7r(p/) 7 (lgg/2) 2 with respect to the 
Z = 28 magic number, that is, a pair of protons align 
into the lgg/2 orbital. This alignment is important for 
the appearance of the shears bands: otherwise no tilted 
rotating state is observed. As for the neutron configura- 
tion, we assume ^(lgg/2) -3 with respect to the N = 50 
magic number. These assumptions lead to a rotational 
band with negative parity, which is consistent with the 
experimental observation. In order to find minima, we 
increase CI = \Cl\ with a step of 0.1 MeV, for each of 
which the minimum in the 9 — <j) plane is searched. Our 
definition of the tilt angles is the same as that in Refs. 
@ and @, that is, ((9,0) = (90°, 0°) corresponds to 
normal collective rotation around the x-axis. This pro- 
cedure is continued until CI = 0.6 MeV. When we come 
to f2 = 0.7 MeV, another minimum appears at which 
the shape is almost spherical. This is consistent with the 
experimental observation of the up-bending seen in 
and B(M1)/B(E2) (see Figs. | and |below), which may 
imply the occurrence of such crossing. Therefore we add 
Cl = 0.65 MeV to the mesh points of the frequency in 
our calculation and the results are shown only up to this 
frequency. Usually we should distinguish the tilt angles 
denned for the rotational frequency vector Cl (denoted 
by 9n and <fin) from those defined for the angular mo- 
mentum vector J (denoted by 9j and </>j). In this study, 
however, we are interested in only each minimum, where 
Cl/ J J and therefore (8q, 4>q) coincides with (9j,4>j), and 
so we simply denote them by (9,4>). The deformations 
/?2 and 7 are calculated from the quadrupole moments 
p2| . We find the deformation slightly decreases from 
(3 2 ~ 0.18 at fi = 0.3 MeV to /3 2 ~ 0.16 at Cl = 0.65 
MeV. These values are close to the fixed value used in 
the P+QQ examination j2C}, that is, e 2 = 0.14. The 
triaxial deformation is at most a few degrees and rather 
small, which can be neglected. This simplifies our calcu- 
lation because we can concentrate on the 2-dimensional 
calculations, where (f> is always set to 0°. Besides, we can 
restrict ourselves to only the range 9 = 0° — 90° thanks 
to the D 2 symmetry |30| , ^l|| . 

Figure [l] shows how the tilt angle 9 at the minima 
changes with respect to the rotational frequency. At low 
rotational frequencies, we find tilted minima appear at 
9 > 50°. With an increase of CI, the tilt angle slightly 
changes but always stays at 9 ~ 55° — 65°. 

In Fig. H the kinematical (J^ 1 **) and dynamical (J^ 2 **) 
moments of inertia are plotted. The experimental values 
are calculated from the transition energies by using the 
finite difference approximations for AI = 1 bands, 
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with Q = 



j( 2 ' does not change so much up to ~ 0.6 MeV. Ex- 
perimental data show an increase of above this fre- 
quency. This might be caused by a crossing as mentioned 
before. Although the calculated is slightly too large, 
as a whole both moments of inertia are well reproduced. 
Small discrepancies seen in Fig. |^ may be cured by in- 
cluding the pairing correlation. 
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FIG. 1. Tilt angle 9 as a function of the rotational fre- 
quency. 
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FIG. 2. Calculated kinematical (J ) and dynamical (J' 2 ') 
moments of inertia together with the experimental data taken 
from Ref. Ed]. Solid and dot-dashed lines represent the cal- 
culated j'~and j' 2 ', respectively. Experimental data are 
shown by filled circles (J' 1 ') and triangles (J' 2 ') connected 
by dotted lines. 



Figure || shows how the direction of the total angular 
momentum vector changes with an increase of the rota- 
tional frequency. Also shown are the net contributions 
from protons and neutrons individually. Because we as- 
sume no core in our calculation, J v and J v are here 
defined as the contributions from all particles below the 
Fermi level, 



Jit — ^ 



N 



J tot J 7T ~\~ J V J 

|J tot | EE y/l(I+l). 

We find, however, almost all contributions come from 
2 proton particles and 3 neutron holes in the lgg/2 or- 
bital. Several proton particles in the lower (p/) orbitals 
also give substantial contributions to the proton angu- 
lar momentum. Because of this, J ^ has not only large 
J^-component but also substantial J z -component even at 
lower frequencies. As the frequency increases, we can see 
that the shears angles decrease with keeping the direc- 
tion of the total angular momentum stayed with nearly 
constant tilt angle. This is explicitly shown in Fig. |^. 
Here the shears angles 0„ and 0t o t are calculated 
from the following semi-classical expressions, 
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FIG. 3. Composition of the total angular momentum at 
Q. = 0.3 MeV (a) and 0.65 MeV (b). J n and J v represent the 
contributions from protons and neutrons, respectively. «7tot 
is the total angular momentum. The angle between «7tot and 
the J z -axis coincides with the tilt angle 8 (see Fig. [l]). 
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Clearly our result shows almost linear decreases of 0^, 
<d u and ©tot- We thus have observed that the shears 
mechanism does appear. 

Another quantity considered as a signal of the shears 
mechanism is the ratio of B(Ml) to B(E2) Q . A smooth 
decrease of the B(M1)/B(E2) ratio is expected in the 
shears bands ||||]. These transition probabilities are 
calculated from the semi-classical expressions For 
-B(Ml), we must calculate the magnetic moments, 
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separately for protons and neutrons. As for the g-factors 
we use the standard values pij] ; gi — 1, g s = 5.58 for 
protons and gi = 0, g s = —3.82 for neutrons, respectively. 
Note that the effective spin g- factor (gi eS ^ ) is equal to the 
free spin g-factor multiplied by 0.7. We do not introduce 
the effective charge for B(E2) because our calculation is 
a fully microscopic one. 
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FIG. 4. The shears angles 8„ and O to t- On is the 
angle between and Jtot, &u between J„ and J tot, and 
Gtot between J ^ and J ' v . 
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FIG. 5. The ratio of B(M1) to B(E2) as a function of the 
rotational frequency. The calculated value is multiplied by 
0.3. Experimental data are taken from Ref. |B0|. 



Figure | shows the B(Ml)/B(E2) ratio. The calculated 
result is attenuated by a factor of 0.3, because only the 
f2 dependence is reliable in the present calculation due to 
the following reasons: 

(1) In our results the effect of the pairing interac- 
tion is not taken into account. As for B(Ml), al- 
most all contributions come from the valence par- 
ticles/holes, while in the case of B(E2) the con- 
tribution from the 'core' part is rather large. Be- 
cause the pairing correlation strongly affects the 
levels near the Fermi level, B(M1) is largely af- 
fected by including the pairing correlation, while 
_B(E2) is not. In our case of 84 Rb, the tendency 
that proton particles in the low-fc levels favor the 
IE-axis (rotational-aligned) and neutron holes in the 
high-Zc levels favor the z-axis (deformation-aligned) 
might be smeared out to some extent by including 
the pairing, because the pairing mixes many levels. 
This can cause the decrease of the shears angle, 
and then B(Ml) would be reduced. Thus we can 
expect that including the pairing will reduce the 
B(M1)/B(E2) ratio. 

(2) The particle-vibration coupling calculations based 
on the RPA show that the in-band S(M1)/B(E2) 
values are reduced if the gamma vibration is col- 
lective enough |35| because it introduces wobbling 
and consequently the overlap between the initial 
and the final states of the transition is reduced. 
This would apply also to the present case, as the 
P+QQ result |2g] suggests that 84 Rb is a gamma- 
soft nucleus. 

The tendency of a smooth decrease of B(M1)/B(E2) up 
to f2 ~ 0.65 MeV is well reproduced as can be seen from 
Fig. ||. Thus our calculation can again show the appear- 
ance of the shears bands. 

To summarize, we have applied the RMF model to 
the newly discovered shears band in 84 Rb. To this 
time only studies based on the P+QQ Hamiltonian had 
been performed as for the mean field approach to the 
shears bands. This is the first examination which is 
based on a more sophisticated mean field model. Our 
RMF calculation shows decreases of the shears angles 
and the B(Ml)/i?(E2) ratio as the frequency increases, 
with keeping nearly constant tilt angle. Thus we could 
reproduce the appearance of the shears mechanism in this 
nucleus. 

Finally, we just give a comment on the case of 82 Rb, 
in which a shears band is observed, too. Unfortunately, 
we could not reproduce the shears band so definitely for 
82 Rb. A possible reason is as follows: In the calculation 
of 82 Rb without pairing, two additional holes must be 
created compared with 84 Rb, that is, there exist 5 neu- 
tron holes with respect to the N = 50 magic number. 
It seems that one of these additional two holes strongly 
favors 8 — 90°, and therefore works in the direction of 



4 



shifting the angular momentum vector toward 6 = 90°. 
Pairing would reduce this effect. 

In the A ~ 80 region, several high-i^ bands have been 
observed |3^ j. It remains unsolved, however, whether 
these are really shears bands or usual high- AT bands in 
the well deformed nuclei. In future, we will examine other 
isotopes in this region, after including the pairing corre- 
lation, to clarify the question above. 
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